Detecting Gene-Gene Interaction among DNA Repair Genes in Chinese non-Syndromic Cleft lip with or Without Palate Trios.

OBJECTIVE: The objective of this study is to investigate the gene-gene interactions associated with NSCL/P among DNA repair genes. DESIGN: This study included 806 NSCL/P case-parent trios from China. Quality control process was conducted for genotyped single nucleotide polymorphisms (SNPs) located in six DNA repair genes (ATR, ERCC4, RFC1, TYMS, XRCC1 and XRCC3). We tested gene-gene interactions with Cordell's method using statistical package TRIO in R software. Bonferroni corrected significance level was set as P = 4.24 × 10-4. We also test the robustness of the interactions by permutation tests. SETTING: Not applicable. PATIENTS/PARTICIPANTS: A total of 806 NSCL/P case-parent trios (complete trios: 682, incomplete trios: 124) with Chinese ancestry. INTERVENTIONS: Not applicable. MAIN OUTCOME MEASURE(S): Not applicable. RESULTS: A total of 118 SNPs were extracted for the interaction tests. Fourteen pairs of significant interactions were identified after Bonferroni correction, which were confirmed in permutation tests. Twelve pairs were between ATR and ERCC4 or XRCC3. The most significant interaction occurred between rs2244500 in TYMS and rs3213403 in XRCC1(P = 8.16 × 10-15). CONCLUSIONS: The current study identified gene-gene interactions among DNA repair genes in 806 Chinese NSCL/P trios, providing additional evidence for the complicated genetic structure underlying NSCL/P. ATR, ERCC4, XRCC3, TYMS and RFC1 were suggested to be possible candidate genes for NSCL/P.

Early Evidence of Chronic Obstructive Pulmonary Disease Obscured by Race-Specific Prediction Equations.

Rationale: The identification of early chronic obstructive pulmonary disease (COPD) is essential to appropriately counsel patients regarding smoking cessation, provide symptomatic treatment, and eventually develop disease-modifying treatments. Disease severity in COPD is defined using race-specific spirometry equations. These may disadvantage non-White individuals in diagnosis and care. Objectives: Determine the impact of race-specific equations on African American (AA) versus non-Hispanic White individuals. Methods: Cross-sectional analyses of the COPDGene (Genetic Epidemiology of Chronic Obstructive Pulmonary Disease) cohort were conducted, comparing non-Hispanic White (n = 6,766) and AA (n = 3,366) participants for COPD manifestations. Measurements and Main Results: Spirometric classifications using race-specific, multiethnic, and "race-reversed" prediction equations (NHANES [National Health and Nutrition Examination Survey] and Global Lung Function Initiative "Other" and "Global") were compared, as were respiratory symptoms, 6-minute-walk distance, computed tomography imaging, respiratory exacerbations, and St. George's Respiratory Questionnaire. Application of different prediction equations to the cohort resulted in different classifications by stage, with NHANES and Global Lung Function Initiative race-specific equations being minimally different, but race-reversed equations moving AA participants to more severe stages and especially between the Global Initiative for Chronic Obstructive Lung Disease (GOLD) stage 0 and preserved ratio impaired spirometry groups. Classification using the established NHANES race-specific equations demonstrated that for each of GOLD stages 1-4, AA participants were younger, had fewer pack-years and more current smoking, but had more exacerbations, shorter 6-minute-walk distance, greater dyspnea, and worse BODE (body mass index, airway obstruction, dyspnea, and exercise capacity) scores and St. George's Respiratory Questionnaire scores. Differences were greatest in GOLD stages 1 and 2. Race-reversed equations reclassified 774 AA participants (43%) from GOLD stage 0 to preserved ratio impaired spirometry. Conclusions: Race-specific equations underestimated disease severity among AA participants. These effects were particularly evident in early disease and may result in late detection of COPD.

Rare variant modifier analysis identifies variants in SEC24D associated with orofacial cleft subtypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p = 6.86 [Formula: see text] 10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p = 0.0009), and created binding sites for 23 TFs, including Pax6 (p = 6.12 [Formula: see text] 10-5) and Pax9 (p = 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study indicates that rare genetic variation may contribute to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Cleft palate (CP) is one of the most common craniofacial birth defects; however, there are relatively few established genetic risk factors associated with its occurrence despite high heritability. Historically, CP has been studied as a single phenotype, although it manifests across a spectrum of defects involving the hard and/or soft palate. We performed a genome-wide association study using transmission disequilibrium tests of 435 case-parent trios to evaluate broad risks for any cleft palate (ACP) (n = 435), and subtype-specific risks for any cleft soft palate (CSP), (n = 259) and any cleft hard palate (CHP) (n = 125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p = 4.24 × 10-8) associated with CHP. One gene at this locus, angiopoietin-like 2 (ANGPTL2), plays a role in osteoblast differentiation. It is expressed both in craniofacial tissue of human embryos and developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p < 5 × 10-6), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios for the 20 loci were most similar across all 3 groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either subtype. We also found nominal evidence of replication (p < 0.05) for 22 SNPs previously associated with orofacial clefts. Our study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

Xanthine oxidoreductase gene polymorphisms are associated with high risk of sepsis and organ failure.

BACKGROUND: Sepsis and associated organ failures confer substantial morbidity and mortality. Xanthine oxidoreductase (XOR) is implicated in the development of tissue oxidative damage in a wide variety of respiratory and cardiovascular disorders including sepsis and sepsis-associated acute respiratory distress syndrome (ARDS). We examined whether single nucleotide polymorphisms (SNPs) in the XDH gene (encoding XOR) might influence susceptibility to and outcome in patients with sepsis. METHODS: We genotyped 28 tag SNPs in XDH gene in the CELEG cohort, including 621 European American (EA) and 353 African American (AA) sepsis patients. Serum XOR activity was measured in a subset of CELEG subjects. Additionally, we assessed the functional effects of XDH variants utilizing empirical data from different integrated software tools and datasets. RESULTS: Among AA patients, six intronic variants (rs206805, rs513311, rs185925, rs561525, rs2163059, rs13387204), in a region enriched with regulatory elements, were associated with risk of sepsis (P < 0.008-0.049). Two out of six SNPs (rs561525 and rs2163059) were associated with risk of sepsis-associated ARDS in an independent validation cohort (GEN-SEP) of 590 sepsis patients of European descent. Two common SNPs (rs1884725 and rs4952085) in tight linkage disequilibrium (LD) provided strong evidence for association with increased levels of serum creatinine (Padjusted<0.0005 and 0.0006, respectively), suggesting a role in increased risk of renal dysfunction. In contrast, among EA ARDS patients, the missense variant rs17011368 (I703V) was associated with enhanced mortality at 60-days (P < 0.038). We found higher serum XOR activity in 143 sepsis patients (54.5 ± 57.1 mU/mL) compared to 31 controls (20.9 ± 12.4 mU/mL, P = 1.96 × 10- 13). XOR activity was associated with the lead variant rs185925 among AA sepsis patients with ARDS (P < 0.005 and Padjusted<0.01). Multifaceted functions of prioritized XDH variants, as suggested by various functional annotation tools, support their potential causality in sepsis. CONCLUSIONS: Our findings suggest that XOR is a novel combined genetic and biochemical marker for risk and outcome in patients with sepsis and ARDS.

Rare variants found in clinical gene panels illuminate the genetic and allelic architecture of orofacial clefting.

PURPOSE: Orofacial clefts (OFCs) are common birth defects including cleft lip, cleft lip and palate, and cleft palate. OFCs have heterogeneous etiologies, complicating clinical diagnostics because it is not always apparent if the cause is Mendelian, environmental, or multifactorial. Sequencing is not currently performed for isolated or sporadic OFCs; therefore, we estimated the diagnostic yield for 418 genes in 841 cases and 294 controls. METHODS: We evaluated 418 genes using genome sequencing and curated variants to assess their pathogenicity using American College of Medical Genetics criteria. RESULTS: 9.04% of cases and 1.02% of controls had "likely pathogenic" variants (P < .0001), which was almost exclusively driven by heterozygous variants in autosomal genes. Cleft palate (17.6%) and cleft lip and palate (9.09%) cases had the highest yield, whereas cleft lip cases had a 2.80% yield. Out of 39 genes with likely pathogenic variants, 9 genes, including CTNND1 and IRF6, accounted for more than half of the yield (4.64% of cases). Most variants (61.8%) were "variants of uncertain significance", occurring more frequently in cases (P = .004), but no individual gene showed a significant excess of variants of uncertain significance. CONCLUSION: These results underscore the etiological heterogeneity of OFCs and suggest sequencing could reduce the diagnostic gap in OFCs.

Causes of and Clinical Features Associated with Death in Tobacco Cigarette Users by Lung Function Impairment.

Rationale: Cigarette smoking contributes to the risk of death through different mechanisms. Objectives: To determine how causes of and clinical features associated with death vary in tobacco cigarette users by lung function impairment. Methods: We stratified current and former tobacco cigarette users enrolled in Genetic Epidemiology of Chronic Obstructive Pulmonary Disease (COPDGene) into normal spirometry, PRISm (Preserved Ratio Impaired Spirometry), Global Initiative for Chronic Obstructive Lung Disease (GOLD) 1-2 COPD, and GOLD 3-4 COPD. Deaths were identified via longitudinal follow-up and Social Security Death Index search. Causes of death were adjudicated after a review of death certificates, medical records, and next-of-kin interviews. We tested associations between baseline clinical variables and all-cause mortality using multivariable Cox proportional hazards models. Measurements and Main Results: Over a 10.1-year median follow-up, 2,200 deaths occurred among 10,132 participants (age 59.5 ± 9.0 yr; 46.6% women). Death from cardiovascular disease was most frequent in PRISm (31% of deaths). Lung cancer deaths were most frequent in GOLD 1-2 (18% of deaths vs. 9-11% in other groups). Respiratory deaths outpaced competing causes of death in GOLD 3-4, particularly when BODE index ⩾7. St. George's Respiratory Questionnaire score ⩾25 was associated with higher mortality in all groups: Hazard ratio (HR), 1.48 (1.20-1.84) normal spirometry; HR, 1.40 (1.05-1.87) PRISm; HR, 1.80 (1.49-2.17) GOLD 1-2; HR, 1.65 (1.26-2.17) GOLD 3-4. History of respiratory exacerbations was associated with higher mortality in GOLD 1-2 and GOLD 3-4, quantitative emphysema in GOLD 1-2, and airway wall thickness in PRISm and GOLD 3-4. Conclusions: Leading causes of death vary by lung function impairment in tobacco cigarette users. Worse respiratory-related quality of life is associated with all-cause mortality regardless of lung function.

Use of the Spirometric "Fixed-Ratio" Underdiagnoses COPD in African-Americans in a Longitudinal Cohort Study.

BACKGROUND: COPD diagnosis is tightly linked to the fixed-ratio spirometry criteria of FEV1/FVC < 0.7. African-Americans are less often diagnosed with COPD. OBJECTIVE: Compare COPD diagnosis by fixed-ratio with findings and outcomes by race. DESIGN: Genetic Epidemiology of COPD (COPDGene) (2007-present), cross-sectional comparing non-Hispanic white (NHW) and African-American (AA) participants for COPD diagnosis, manifestations, and outcomes. SETTING: Multicenter, longitudinal US cohort study. PARTICIPANTS: Current or former smokers with ≥ 10-pack-year smoking history enrolled at 21 clinical centers including over-sampling of participants with known COPD and AA. Exclusions were pre-existing non-COPD lung disease, except for a history of asthma. MEASUREMENTS: Subject diagnosis by conventional criteria. Mortality, imaging, respiratory symptoms, function, and socioeconomic characteristics, including area deprivation index (ADI). Matched analysis (age, sex, and smoking status) of AA vs. NHW within participants without diagnosed COPD (GOLD 0; FEV1 ≥ 80% predicted and FEV1/FVC ≥ 0.7). RESULTS: Using the fixed ratio, 70% of AA (n = 3366) were classified as non-COPD, versus 49% of NHW (n = 6766). AA smokers were younger (55 vs. 62 years), more often current smoking (80% vs. 39%), with fewer pack-years but similar 12-year mortality. Density distribution plots for FEV1 and FVC raw spirometry values showed disproportionate reductions in FVC relative to FEV1 in AA that systematically led to higher ratios. The matched analysis demonstrated GOLD 0 AA had greater symptoms, worse DLCO, spirometry, BODE scores (1.03 vs 0.54, p < 0.0001), and greater deprivation than NHW. LIMITATIONS: Lack of an alternative diagnostic metric for comparison. CONCLUSIONS: The fixed-ratio spirometric criteria for COPD underdiagnosed potential COPD in AA participants when compared to broader diagnostic criteria. Disproportionate reductions in FVC relative to FEV1 leading to higher FEV1/FVC were identified in these participants and associated with deprivation. Broader diagnostic criteria for COPD are needed to identify the disease across all populations.

Rare genetic variants in SEC24D modify orofacial cleft phenotypes.

As one of the most common structural birth defects, orofacial clefts (OFCs) have been studied for decades, and recent studies have demonstrated that there are genetic differences between the different phenotypic presentations of OFCs. However, the contribution of rare genetic variation genome-wide to different subtypes of OFCs has been understudied, with most studies focusing on common genetic variation or rare variation within targeted regions of the genome. Therefore, we used whole-genome sequencing data from the Gabriella Miller Kids First Pediatric Research Program to conduct a gene-based burden analysis to test for genetic modifiers of cleft lip (CL) vs cleft lip and palate (CLP). We found that there was a significantly increased burden of rare variants in SEC24D in CL cases compared to CLP cases (p=6.86×10-7). Of the 15 variants within SEC24D, 53.3% were synonymous, but overlapped a known craniofacial enhancer. We then tested whether these variants could alter predicted transcription factor binding sites (TFBS), and found that the rare alleles destroyed binding sites for 9 transcription factors (TFs), including Pax1 (p=0.0009), and created binding sites for 23 TFs, including Pax6 (p=6.12×10-5) and Pax9 (p= 0.0001), which are known to be involved in normal craniofacial development, suggesting a potential mechanism by which these synonymous variants could have a functional impact. Overall, this study demonstrates that rare genetic variation contributes to the phenotypic heterogeneity of OFCs and suggests that regulatory variation may also contribute and warrant further investigation in future studies of genetic variants controlling risk to OFC.

Trio-based GWAS identifies novel associations and subtype-specific risk factors for cleft palate.

Orofacial clefts (OFCs) are the most common craniofacial birth defects and are often categorized into two etiologically distinct groups: cleft lip with or without cleft palate (CL/P) and isolated cleft palate (CP). CP is highly heritable, but there are still relatively few established genetic risk factors associated with its occurrence compared to CL/P. Historically, CP has been studied as a single phenotype despite manifesting across a spectrum of defects involving the hard and/or soft palate. We performed GWAS using transmission disequilibrium tests using 435 case-parent trios to evaluate broad risks for any cleft palate (ACP, n=435), as well as subtype-specific risks for any cleft soft palate (CSP, n=259) and any cleft hard palate (CHP, n=125). We identified a single genome-wide significant locus at 9q33.3 (lead SNP rs7035976, p=4.24×10 -8 ) associated with CHP. One gene at this locus, angiopoietin-like 2 ( ANGPTL2 ), plays a role in osteoblast differentiation. It is expressed in craniofacial tissue of human embryos, as well as in the developing mouse palatal shelves. We found 19 additional loci reaching suggestive significance (p<5×10 -6 ), of which only one overlapped between groups (chromosome 17q24.2, ACP and CSP). Odds ratios (ORs) for each of the 20 loci were most similar across all three groups for SNPs associated with the ACP group, but more distinct when comparing SNPs associated with either the CSP or CHP groups. We also found nominal evidence of replication (p<0.05) for 22 SNPs previously associated with cleft palate (including CL/P). Interestingly, most SNPs associated with CL/P cases were found to convey the opposite effect in those replicated in our dataset for CP only. Ours is the first study to evaluate CP risks in the context of its subtypes and we provide newly reported associations affecting the broad risk for CP as well as evidence of subtype-specific risks.

X chromosome associations with chronic obstructive pulmonary disease and related phenotypes: an X chromosome-wide association study.

BACKGROUND: The association between genetic variants on the X chromosome to risk of COPD has not been fully explored. We hypothesize that the X chromosome harbors variants important in determining risk of COPD related phenotypes and may drive sex differences in COPD manifestations. METHODS: Using X chromosome data from three COPD-enriched cohorts of adult smokers, we performed X chromosome specific quality control, imputation, and testing for association with COPD case-control status, lung function, and quantitative emphysema. Analyses were performed among all subjects, then stratified by sex, and subsequently combined in meta-analyses. RESULTS: Among 10,193 subjects of non-Hispanic white or European ancestry, a variant near TMSB4X, rs5979771, reached genome-wide significance for association with lung function measured by FEV1/FVC ([Formula: see text] 0.020, SE 0.004, p 4.97 × 10-08), with suggestive evidence of association with FEV1 ([Formula: see text] 0.092, SE 0.018, p 3.40 × 10-07). Sex-stratified analyses revealed X chromosome variants that were differentially trending in one sex, with significantly different effect sizes or directions. CONCLUSIONS: This investigation identified loci influencing lung function, COPD, and emphysema in a comprehensive genetic association meta-analysis of X chromosome genetic markers from multiple COPD-related datasets. Sex differences play an important role in the pathobiology of complex lung disease, including X chromosome variants that demonstrate differential effects by sex and variants that may be relevant through escape from X chromosome inactivation. Comprehensive interrogation of the X chromosome to better understand genetic control of COPD and lung function is important to further understanding of disease pathology. Trial registration Genetic Epidemiology of COPD Study (COPDGene) is registered at ClinicalTrials.gov, NCT00608764 (Active since January 28, 2008). Evaluation of COPD Longitudinally to Identify Predictive Surrogate Endpoints Study (ECLIPSE), GlaxoSmithKline study code SCO104960, is registered at ClinicalTrials.gov, NCT00292552 (Active since February 16, 2006). Genetics of COPD in Norway Study (GenKOLS) holds GlaxoSmithKline study code RES11080, Genetics of Chronic Obstructive Lung Disease.

Evidence of the folate-mediated one-carbon metabolism pathway genes in controlling the non-syndromic oral clefts risks.

The folate-mediated one-carbon metabolism pathway is thought to play an important role in the etiology of non-syndromic oral clefts (NSOFC), although none of the genes in this pathway has shown significant signals in genome-wide association studies (GWAS). Recent evidence indicated that enhanced understanding could be gained by aggregating multiple SNPs effect simultaneously into polygenic risk score (PRS) to assess its association with disease risks. This study is aimed to assess the association between the genetic effect of folate-mediated one-carbon metabolism pathway and NSOFC risks using PRS based on a case-parent trio design. A total of 297 SNPs mapped from 18 genes in the folate-mediated one-carbon metabolism pathway were aggregated from a GWAS of 2458 case-parent trios recruited from an international consortium. We found a PRS based on the folate-mediated one-carbon metabolism pathway was significant among all NSOFC trios (OR = 1.95, 95% CI: 1.66-2.28, p = 2.39 × 10-16 ), as well as two major subtypes, non-syndromic cleft lip with or without cleft palate (NSCL/P) trios (OR = 1.71, 95% CI: 1.50-1.96, p = 7.66 × 10-15 ) and non-syndromic cleft palate only (NSCPO) trios (OR = 1.51, 95% CI: 1.36-1.68, p = 2.1 × 10-14 ). Similar results were also observed in further subgroup analyses stratified into Asian and European trios. The averaged PRS of the folate-mediated one-carbon metabolism pathway varied between the NSOFC case group and its comparison group (p < 0.05) with higher average PRS in the cases. Moreover, the top 5% pathway PRS group had 2.25 (95% CI: 1.85-2.73) times increased NSOFC risk, also 3.09 (95% CI: 2.50-3.81) and 2.06 (95% CI: 1.39-3.02) times increased risk of NSCL/P and NSCPO compared to the remainder of the distribution. The results of our study confirmed the folate-mediated one-carbon metabolism pathway was important in controlling risk to NSOFC and this study enhanced evidence towards understanding the genetic risks of NSOFC.

Damaging Mutations in AFDN Contribute to Risk of Nonsyndromic Cleft Lip With or Without Cleft Palate.

Novel or rare damaging mutations have been implicated in the developmental pathogenesis of nonsyndromic cleft lip with or without cleft palate (nsCL ± P). Thus, we investigated the human genome for high-impact mutations that could explain the risk of nsCL ± P in our cohorts.We conducted next-generation sequencing (NGS) analysis of 130 nsCL ± P case-parent African trios to identify pathogenic variants that contribute to the risk of clefting. We replicated this analysis using whole-exome sequence data from a Brazilian nsCL ± P cohort. Computational analyses were then used to predict the mechanism by which these variants could result in increased risks for nsCL ± P.We discovered damaging mutations within the AFDN gene, a cell adhesion molecule (CAMs) that was previously shown to contribute to cleft palate in mice. These mutations include p.Met1164Ile, p.Thr453Asn, p.Pro1638Ala, p.Arg669Gln, p.Ala1717Val, and p.Arg1596His. We also discovered a novel splicing p.Leu1588Leu mutation in this protein. Computational analysis suggests that these amino acid changes affect the interactions with other cleft-associated genes including nectins (PVRL1, PVRL2, PVRL3, and PVRL4) CDH1, CTNNA1, and CTNND1.This is the first report on the contribution of AFDN to the risk for nsCL ± P in humans. AFDN encodes AFADIN, an important CAM that forms calcium-independent complexes with nectins 1 and 4 (encoded by the genes PVRL1 and PVRL4). This discovery shows the power of NGS analysis of multiethnic cleft samples in combination with a computational approach in the understanding of the pathogenesis of nsCL ± P.

Whole-genome sequencing reveals de-novo mutations associated with nonsyndromic cleft lip/palate.

The majority (85%) of nonsyndromic cleft lip with or without cleft palate (nsCL/P) cases occur sporadically, suggesting a role for de novo mutations (DNMs) in the etiology of nsCL/P. To identify high impact protein-altering DNMs that contribute to the risk of nsCL/P, we conducted whole-genome sequencing (WGS) analyses in 130 African case-parent trios (affected probands and unaffected parents). We identified 162 high confidence protein-altering DNMs some of which are based on available evidence, contribute to the risk of nsCL/P. These include novel protein-truncating DNMs in the ACTL6A, ARHGAP10, MINK1, TMEM5 and TTN genes; as well as missense variants in ACAN, DHRS3, DLX6, EPHB2, FKBP10, KMT2D, RECQL4, SEMA3C, SEMA4D, SHH, TP63, and TULP4. Many of these protein-altering DNMs were predicted to be pathogenic. Analysis using mouse transcriptomics data showed that some of these genes are expressed during the development of primary and secondary palate. Gene-set enrichment analysis of the protein-altering DNMs identified palatal development and neural crest migration among the few processes that were significantly enriched. These processes are directly involved in the etiopathogenesis of clefting. The analysis of the coding sequence in the WGS data provides more evidence of the opportunity for novel findings in the African genome.

A genome-wide association study of bronchodilator response in participants of European and African ancestry from six independent cohorts.

Introduction: Bronchodilator response (BDR) is a measurement of acute bronchodilation in response to short-acting ß2-agonists, with a heritability between 10 and 40%. Identifying genetic variants associated with BDR may lead to a better understanding of its complex pathophysiology. Methods: We performed a genome-wide association study (GWAS) of BDR in six adult cohorts with participants of European ancestry (EA) and African ancestry (AA) including community cohorts and cohorts ascertained on the basis of obstructive pulmonary disease. Validation analysis was carried out in two paediatric asthma cohorts. Results: A total of 10 623 EA and 3597 AA participants were included in the analyses. No single nucleotide polymorphism (SNP) was associated with BDR at the conventional genome-wide significance threshold (p<5×10-8). Performing fine mapping and using a threshold of p<5×10-6 to identify suggestive variants of interest, we identified three SNPs with possible biological relevance: rs35870000 (within FREM1), which may be involved in IgE- and IL5-induced changes in airway smooth muscle cell responsiveness; rs10426116 (within ZNF284), a zinc finger protein, which has been implicated in asthma and BDR previously; and rs4782614 (near ATP2C2), involved in calcium transmembrane transport. Validation in paediatric cohorts yielded no significant SNPs, possibly due to age-genotype interaction effects. Conclusion: Ancestry-stratified and ancestry-combined GWAS meta-analyses of over 14 000 participants did not identify genetic variants associated with BDR at the genome-wide significance threshold, although a less stringent threshold identified three variants showing suggestive evidence of association. A common definition and protocol for measuring BDR in research may improve future efforts to identify variants associated with BDR.

Benchmarking statistical methods for analyzing parent-child dyads in genetic association studies.

Genetic association studies of child health outcomes often employ family-based study designs. One of the most popular family-based designs is the case-parent trio design that considers the smallest possible nuclear family consisting of two parents and their affected child. This trio design is particularly advantageous for studying relatively rare disorders because it is less prone to type 1 error inflation due to population stratification compared to population-based study designs (e.g., case-control studies). However, obtaining genetic data from both parents is difficult, from a practical perspective, and many large studies predominantly measure genetic variants in mother-child dyads. While some statistical methods for analyzing parent-child dyad data (most commonly involving mother-child pairs) exist, it is not clear if they provide the same advantage as trio methods in protecting against population stratification, or if a specific dyad design (e.g., case-mother dyads vs. case-mother/control-mother dyads) is more advantageous. In this article, we review existing statistical methods for analyzing genome-wide marker data on dyads and perform extensive simulation experiments to benchmark their type I errors and statistical power under different scenarios. We extend our evaluation to existing methods for analyzing a combination of case-parent trios and dyads together. We apply these methods on genotyped and imputed data from multiethnic mother-child pairs only, case-parent trios only or combinations of both dyads and trios from the Gene, Environment Association Studies consortium (GENEVA), where each family was ascertained through a child affected by nonsyndromic cleft lip with or without cleft palate. Results from the GENEVA study corroborate the findings from our simulation experiments. Finally, we provide recommendations for using statistical genetic association methods for dyads.

Individual and Combined Association Between Prenatal Polysubstance Exposure and Childhood Risk of Attention-Deficit/Hyperactivity Disorder.

Importance: Polysubstance use among pregnant women has increased because of the opioid epidemic and the increasing legalization of cannabis along with persistent tobacco and alcohol consumption. Previous research on prenatal substance use and the child's risk of attention-deficit/hyperactivity disorder (ADHD) has mostly focused on single-substance exposures; simultaneous examination of multiple substance use and assessment of their synergistic health consequences is needed. Objectives: To assess the consequences of the use of specific substances during pregnancy, investigate whether the interaction of multiple prenatal substance exposures is associated with increases in the risk of childhood ADHD, and estimate the aggregate burden of polysubstance exposure during gestation. Design, Setting, and Participants: This cohort study analyzed data from the Boston Birth Cohort from 1998 to 2019. The sample of the present study comprised a multiethnic urban cohort of mother-child pairs who were predominantly low income. A total of 3138 children who were enrolled shortly after birth at Boston Medical Center were included and followed up from age 6 months to 21 years. Exposures: Substance use during pregnancy was identified based on self-reported tobacco smoking, alcohol consumption, and use of cannabis, cocaine, or opioids in any trimester of pregnancy. Diagnostic codes for neonatal opioid withdrawal syndrome or neonatal abstinence syndrome from the International Classification of Diseases, Ninth Revision, and the International Classification of Diseases, Tenth Revision, were also used to identify opioid exposure during gestation. Main Outcomes and Measures: ADHD diagnosis in the child's electronic medical record. Results: Among 3138 children (1583 boys [50.4%]; median age, 12 years [IQR, 9-14 years]; median follow-up, 10 years [IQR, 7-12 years]) in the final analytic sample, 486 (15.5%) had an ADHD diagnosis and 2652 (84.5%) were neurotypical. The median postnatal follow-up duration was 12 years (IQR, 9-14 years). Among mothers, 46 women (1.5%) self-identified as Asian (non-Pacific Islander), 701 (22.3%) as Hispanic, 1838 (58.6%) as non-Hispanic Black, 227 (7.2%) as non-Hispanic White, and 326 (10.4%) as other races and/or ethnicities (including American Indian or Indigenous, Cape Verdean, Pacific Islander, multiracial, other, or unknown). A total of 759 women (24.2%) reported the use of at least 1 substance during pregnancy, with tobacco being the most frequently reported (580 women [18.5%]). Cox proportional hazards models revealed that opioid exposure (60 children) had the highest adjusted hazard ratio (HR) for ADHD (2.19; 95% CI, 1.10-4.37). After including main statistical effects of all individual substances in an elastic net regression model, the HR of opioids was reduced to 1.60, and evidence of a statistical interaction between opioids and both cannabis and alcohol was found, producing 1.42 and 1.15 times higher risk of ADHD, respectively. The interaction between opioids and smoking was also associated with a higher risk of ADHD (HR, 1.17). Conclusions and Relevance: The findings of this study suggest that it is important to consider prenatal concurrent exposure to multiple substances and their possible interactions when counseling women regarding substance use during pregnancy, the future risk of ADHD for their children, and strategies for cessation and treatment programs.

Accurate assignment of disease liability to genetic variants using only population data.

PURPOSE: The growing size of public variant repositories prompted us to test the accuracy of pathogenicity prediction of DNA variants using population data alone. METHODS: Under the a priori assumption that the ratio of the prevalence of variants in healthy population vs that in affected populations form 2 distinct distributions (pathogenic and benign), we used a Bayesian method to assign probability to a variant belonging to either distribution. RESULTS: The approach, termed Bayesian prevalence ratio (BayPR), accurately parsed 300 of 313 expertly curated CFTR variants: 284 of 296 pathogenic/likely pathogenic variants in 1 distribution and 16 of 17 benign/likely benign variants in another. BayPR produced an area under the receiver operating characteristic curve of 0.99 for 103 functionally confirmed missense CFTR variants, which is equal to or exceeds 10 commonly used algorithms (area under the receiver operating characteristic curve range = 0.54-0.99). Application of BayPR to expertly curated variants in 8 genes associated with 7 Mendelian conditions led to the assignment of a disease-causing probability of ≥80% to 1350 of 1374 (98.3%) pathogenic/likely pathogenic variants and of ≤20% to 22 of 23 (95.7%) benign/likely benign variants. CONCLUSION: Irrespective of the variant type or functional effect, the BayPR approach provides probabilities of pathogenicity for DNA variants responsible for Mendelian disorders using only the variant counts in affected and unaffected population samples.

Whole genome sequence analysis of platelet traits in the NHLBI Trans-Omics for Precision Medicine (TOPMed) initiative.

Platelets play a key role in thrombosis and hemostasis. Platelet count (PLT) and mean platelet volume (MPV) are highly heritable quantitative traits, with hundreds of genetic signals previously identified, mostly in European ancestry populations. We here utilize whole genome sequencing (WGS) from NHLBI's Trans-Omics for Precision Medicine initiative (TOPMed) in a large multi-ethnic sample to further explore common and rare variation contributing to PLT (n = 61 200) and MPV (n = 23 485). We identified and replicated secondary signals at MPL (rs532784633) and PECAM1 (rs73345162), both more common in African ancestry populations. We also observed rare variation in Mendelian platelet-related disorder genes influencing variation in platelet traits in TOPMed cohorts (not enriched for blood disorders). For example, association of GP9 with lower PLT and higher MPV was partly driven by a pathogenic Bernard-Soulier syndrome variant (rs5030764, p.Asn61Ser), and the signals at TUBB1 and CD36 were partly driven by loss of function variants not annotated as pathogenic in ClinVar (rs199948010 and rs571975065). However, residual signal remained for these gene-based signals after adjusting for lead variants, suggesting that additional variants in Mendelian genes with impacts in general population cohorts remain to be identified. Gene-based signals were also identified at several genome-wide association study identified loci for genes not annotated for Mendelian platelet disorders (PTPRH, TET2, CHEK2), with somatic variation driving the result at TET2. These results highlight the value of WGS in populations of diverse genetic ancestry to identify novel regulatory and coding signals, even for well-studied traits like platelet traits.